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SECTION  1 
INTRODUCTION 


1.1  BACKGROUND. 

There  has  been  renewed  interest  in  explosions  inside  chambers  with  vents  for  both  safety 
and  military  reasons.  Escape  of  hazardous  vented  materials  is  of  great  concern  in  ensuring 
civilian  safety;  it  is  also  quite  important  in  estimating  collateral  effects  from  military  opera¬ 
tions  against  stored  chemical/biological  agents.  The  effects  of  internal  blast  on  the  venting 
of  material  contained  within  is  currently  not  fully  understood.  As  such,  more  experiments 
are  being  conducted  to  understand  venting  from  internal  blast  and  new  predictive  models 
are  being  developed.  In  designing  blast-venting  experiments,  there  appears  to  be  a  need 
for  a  simple  model  for  rapid  trade-off  studies  that  combines  the  effects  of  fluid  mechanics, 
thermal  and  mechanical  factors  on  venting.  In  this  report,  a  blast  venting  model  for  a 
single  chamber  is  described. 

1.2  SCOPE  OF  THE  REPORT. 

This  report  documents  a  venting  model  embodied  in  a  computer  code  named  MATVENT, 
and  several  calculations  carried  out  using  the  code  to  learn  about  the  relative  importance 
of  various  factors  which  could  affect  venting;  of  particular  interest  are  the  duration  of 
venting  and  the  total  mass  vented.  The  model  considers  the  venting  of  gases  from  a  blast 
inside  a  single,  mechanically  and  thermally  responding  chamber  with  a  small  opening 
to  the  ambient  atmosphei'e  (Figure  1-1).  “Mechanically  Responding”  means  that  the 
interior  walls  of  the  chamber  I'espond  to  chamber  pressure.  “Thermally  responding”  refers 
to  the  change  in  the  energetics  and  mechanics  of  venting  because  of  the  loss  of  blast 
energy  through  the  walls  by  thermal  effects  such  as  conduction  and  convection;  thermal 
response  also  includes  the  volume  changes  brought  about  because  of  thermal  expansion  or 
contraction  of  the  walls. 

The  essentials  of  the  MATVENT  venting  model  are  detailed  in  Section  2,  along  with  a 
list  of  the  needed  input  data  and  other  assumptions.  Section  3  gives  a  brief  description 
of  the  code  and  the  high-level  language  MATLAB  in  which  the  code  is  written.  Section  4 
describes  the  results  from  several  parametric  MATVENT  calculations  employing  the  same 
geometry  and  explosion  yield  as  a  recent  experiment,  briefly  described  in  Section  4.1.  For 
the  same  explosive  yield,  chamber  contents  and  chamber/vent  dimensions,  the  effects  on 
venting  of: 
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Figure  1-1.  Schematic  for  the  blast  venting  problem. 

•  assuming  different  specific  heat  ratios  7 

•  assuming  diffei'ent  fluidized  mass  fractions  of  the  chamber  contents  in  a  rigid  and 
insulated  chamber.  The  term  “fluidized”  refers  to,  the  “dense  -gas”  analytic  approx¬ 
imation  which  treats  the  mass  of  loose  solid  material  as  if  it  were  a  gas  uniformly 
distributed  over  the  volume,  intimately  mixed  and  in  equilibrium  with  the  chamber 
air  and  explosion  gases. 

•  assuming  a  chamber  which  is  insulated  but  deformable 

•  assuming  a  chamber  which  loses  energy  by  thermal  conduction  only,  but  does  not 
respond  to  pressure 

•  assuming  a  chamber  which  loses  energy  by  convection  and  conduction  effects,  without 
responding  to  pressure,  and 

•  assuming  a  chamber  which  is  fully  responsive  to  both  pressure  and  thermal  effects. 
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are  investigated  parametrically  with  the  objective  of  recognizing  the  relative  importance  of 
the  assumptions  in  venting  estimates.  Additional  calculations  to  resolve  the  applicability 
of  yield-scaling,  important  in  the  design  of  scaled  experiments,  are  reported.  The  effect  of 
each  set  of  assumptions  on  the  important  quantities,  such  as  venting  rates,  total  vented 
mass,  and  total  venting  time  are  obtained  numerically.  Section  5  discusses  the  calculational 
results  and  the  limitations.  Section  6  summarizes  the  results  and  provides  recommenda¬ 
tions  for  further  work  to  improve  the  model.  Appendix  A  provides  the  model  equations 
which  have  been  coded  in  MATVENT.  In  Appendbc  B,  a  listing  of  the  basic  MATVENT 
code  is  provided  along  with  the  results  of  a  run  from  an  example  case.  Finally,  Appendix 
C  provides  definitions  for  the  symbols  used  in  the  body  of  the  report. 
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SECTION  2 

DESCRIPTION  OF  THE  MODEL 

This  section  provides  descriptions  of  the  elements  of  the  basic  version  of  the  MATVENT 
model  as  well  as  modifications  for  parametric  studies.  The  input  data  needs,  the  fluid  me¬ 
chanics  model  element,  the  mechajiical  response  element  and  the  thermal  response  element 
are  detailed. 

2.1  INPUT  DATA  NEEDS  AND  OTHER  ASSUMPTIONS. 

The  parameters  needed  as  inputs  to  the  model  and  assumed  as  given,  are  the  following: 

(a)  The  thermal  and  mechanical  properties  of  the  material  out  of  which  the  chamber  is 
built  (steel,  concrete,  etc.); 

(b)  The  convective  heat  transfer  coefficient  ‘u’  (also  referred  to  sometimes  as  “unit 
surface  conductance”); 

(c)  The  volume  and  linear  dimensions  of  the  chamber; 

(d)  The  cross-sectional  area  of  the  vent; 

(e)  The  mass  of  loose  material  contained  within  the  chamber; 

(f)  The  mass  of  explosive,  and  the  effective  yield  of  the  explosive  per  unit  mass; 

(g)  The  fraction  /  of  loose  material  fluidized  by  the  explosion;  and 

(h)  The  specific  heats  and  specific  heat  ratio  (7)  of  the  chamber  gases  subsequent  to  the 
explosion. 

2.2  MODEL  FOR  THE  FLUID  MECHANICS  OF  VENTING. 

The  model  assumes  the  fluid  mechanical  process  inside  the  chamber  to  be  quasi-static 
and  governed  by  the  adiabatic  relations  for  a  calorically  and  thermally  perfect  (“dense”) 
gas,  as  is  usually  done  (Reference  1)  for  flow  out  of  a  reservoir  with  a  vent.  The  model 
does  not  treat  the  initial  shock  reverberation  period  which  might  last  for  a  very  short 
time  (of  the  order  of  tens  of  milliseconds)  compared  to  total  venting  durations.  Time  zero 
in  MATVENT  can  be  shifted  if  need  be  to  a  time  marking  the  end  of  the  initial  shock 
reverberation  intex'val.  The  shock  reverberation  interval  is  typically  a  few  shock  trausit 
times  across  a  representative  linear  dimension  of  the  chamber,  such  as  the  hydraulic  mean 
diameter. 
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Isentropic  relations  applicable  under  the  quasi-static  assumption  enable  the  calculation 
of  the  time-dependent  chamber  pressure,  temperature  and  density,  as  well  as  the  vented 
fluidized  mass  and  venting  duration,  in  terms  of  the  initial  conditions  obtained  by  dis¬ 
tributing  the  energy  of  the  explosion  and  the  mass  of  fluidized  material  uniformly  over  the 
volume  of  the  chamber.  The  equations  and  other  details  of  fluid  mechanical  modeling  are 
given  in  Appendix  A. 

2.3  CHAMBER  RESPONSE. 

The  modeling  for  mechanical  and  thermal  response  assumes  that  the  time  scale  for  fluid 
mechanical  equilibration  is  much  smaller  than  those  for  mechanical  or  thermal  response. 
Even  where  this  assumption  is  not  quite  valid,  the  model  will  still  provide  a  “feel”  for  the 
effects  of  such  response  on  the  venting,  and  thus  help  orient  experiment  design. 

2.3.1  Mechanical  Response  Model. 

The  mechanical  response  of  the  chamber  to  the  explosion  inside  is  modeled  in  the  following 
simple  manner.  At  each  timestep  6t  of  the  calculation,  the  incremental  overpressure  im¬ 
pulse  p-6t  is  obtained  and  divided  by  the  specific  acoustic  impedance  w  of  the  wall  material 
to  obtain  an  incremental  displacement.  The  volume  change  increment  SV  is  the  product  of 
the  total  wall  area  ARq  and  the  incremental  displacement.  The  total  volume  change  at  any 
time  is  the  sum  of  the  incremental  volume  changes  up  to  that  time.  Thennal  expansion 
can,  as  described  in  Section  2.3.2,  also  be  included  at  each  timestep,  leading  to  a  ther¬ 
mal  modification  (reduction)  of  the  purely  mechanical  inci'emental  volume  change.  The 
updated  volume  at  any  time  is  obtained  by  adding  the  volume  change  during  a  timestep 
to  the  volume  at  the  beginning  of  the  timestep  ( “old  volume” .  The  density  and  pressure 
inside  the  chamber  are  adjusted  using  the  ratio  of  the  updated  and  “old”  volumes.  The 
temperatui'e  ratio  is  also  adjusted  to  an  updated  value.  The  equations  describing  modeling 
for  mechanical  response  are  detailed  in  Appendbc  A. 

2.3.2  Thermal  Response  Models. 

Thermal  response  was  initially  modeled  as  if  conduction  alone  were  important.  Later  a 
model,  presumed  more  realistic,  was  developed  to  teike  both  convection  near  the  walls  as 
well  as  conduction  through  the  walls  into  account. 
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2.3.2. 1  Thermal  Response  Model  For  Conduction  Alone.  In  modeling  the  thermal  re¬ 
sponse  alone,  first  a  model  for  heat  loss  by  conduction  only  was  formulated.  In  this 
“conduction  only”  model,  an  “effective”  penetration  thickness  ‘^a  •  t'  for  the  heat  pulse 
was  calculated  based  on  the  thermal  diffusivity  a  of  the  chamber  material  and  the  time 
after  blast  t.  This  “effective  thickness”  is  a  fictitious  thickness  of  the  wall  which,  for  the 
purpose  of  heat  loss  calculation  only,  is  assumed  to  rise  instantaneously  to  the  applied 
(time-dependent)  chamber  temperature.  The  total  heat  loss  Q  to  the  chamber  walls  at 
time  t  is  the  product  of  the  penetration  thickness,  the  wall  surface  area,  the  density,  the 
heat  capacity  and  the  chamber  temperature  in  excess  of  ambient: 


Q  =  ~Pw  •  •  y/a.-  Vi  •  C  •  [Tg  —  Toa) 


(2.1) 


By  differentiating  the  expression  for  heat  loss  at  any  time,  the  incremental  heat  loss  6Q 
over  the  timestep  6t  can  be  obtained: 


6Q  —  — pvv  •  ARo  y/Si-C 


'^Oa)  +  VtTg 


■6t 


(2.2) 


The  heat  loss  increment  6Q  can  be  expressed  in  terms  of  the  chamber  pressure  ratio  Xg{t) 
and  its  time  derivative  Xg{t): 


6Q  —  pyf  •  ARq  •  y^a,  •  s/t  •  C  •  Too  • 

The  net  heat  loss  up  to  any  time  is  the  sum  of  the  incremental  losses  over  the  timesteps;  the 
net  heat  loss  divided  by  the  mass  of  the  fluidized  material  in  the  chamber  and  the  specific 
heat  of  the  “dense”  gas  results  in  the  change  (decrease)  of  the  chamber  temperature  due 
to  conduction  only.  The  modified  chamber  temperature  is  then  used  to  calculate  modified 
pressure  and  density,  now  adjusted  for  heat  loss.  The  calculation  proceeds  until  there 
is  pressure  equilibration  between  the  chamber  and  the  ambient  air  and  there  is  no  more 
efflux.  The  chamber  temperature  generally  remains  higher  than  the  ambient  at  this  point, 
while  the  density  is  lower.  Details  of  the  modeling  for  thermal  response  due  to  conduction 
only  are  given  in  Appendix  A.  A  version  of  MATVENT  was  run  for  the  case  of  thermal 
response  due  to  conduction  only  in  order  to  obtain  a  sense  of  how  venting  might  be  affected 
by  conduction  alone.  However,  thermal  response  by  both  convection  and  conduction,  as 
in  Section  2.3. 2. 2  was  deemed  more  realistic  and  is  included  in  the  code  listing,  Appendix 
B,  rather  than  response  by  conduction  only. 


x„  •  Xc 


St  (2.3) 
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2. 3.2.2  Thermal  Response  Model  For  Convection  and  Conduction.  The  temperature  of 
the  wall  surface  in  contact  with  the  chamber  gases  is  in  practicality,  not  the  same  eis  the 
chamber  interior  temperature.  A  drop  in  the  gas  temperature  in  a  thin  layer  adjacent  to 
the  wall,  from  the  interior  chamber  temperature  Tg  to  the  temperature  of  the  wall  surface 
Ty,  exists  in  general,  due  to  convection  in  the  chamber.  A  simple  model  for  heat  loss 
by  both  convection  and  conduction  was  therefore  developed  in  order  to  account  for  both 
the  thermal  resistance  at  the  chamber  walls  and  the  heat  conduction  in  the  body  of  the 
chamber  material.  This  appears  to  be  a  more  realistic  approach  to  modeling.  In  this  model, 
the  “effective”  thickness  of  the  chamber  was  assumed  to  rise  to  a  temperature  lower  than 
the  chamber  temperature  at  any  given  time.  The  temperature  of  the  “effective”  thickness 
was  calculated  from  a  modified  Newtonian  heat  loss  equation  where  the  heat  transfer  rate 
across  the  surface  layer  adjacent  to  the  wall  is  in  balance  with  the  rate  of  (time-dependent) 
heat  gain  by  the  effective  thickness  (Reference  2): 

g  =  Pw-ARo  =  u- ARo-  (Tp-Tw)  (2.4) 

Solving  the  above  differential  equation  with  the  initial  condition  Ty,  —  Toa  at  t  =  0  (and 
noting  that  the  chamber  temperature  Tg  is  to  be  regarded  as  invariant  in  keeping  with  the 
quasi-static  assumption  for  the  chamber  conditions)  leads  to  the  wall  temperature  Tw  at 
time  t  and  the  heat  loss  increment  6Q  over  a  timestep  St  as: 

T^  -  Toa  =  {Tg  -  Toa)  •  (1  -  e-2€'^) 

<5Q  =  -u  •  ARc  {Tg  -  Toa)  •  •  6t  (2.5) 

The  net  heat  loss  up  to  any  time  is  calculated  as  the  sum  of  the  incremental  losses  over 
the  elapsed  timesteps.  The  rest  of  the  procedure  follows  as  in  the  previous  Section  2.3.2. 1. 

The  volume  decrease  produced  in  the  chamber  purely  due  to  thermal  expansion  of  the 
“effective”  thickness  can  be  formulated  as  follows: 

=  —ARq  •  €  •  ■  y/i  •  (Tw  —  7ba)  •  St 

=  -ARo  .e-^a-Vi-{Tg- Tq^)  ■  (l  - e-^^^)  •  St  (2.6) 

The  equations  in  the  modeling  for  thermal  response  by  both  convection  and  conduction 
are  detailed  in  Appendix  A. 
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2.3.3  Combined  Response  Model. 

Both  mechanical  and  thermal  responses  can  simultaneously  be  taken  into  account  in  the 
“full  response”  model  as  given  in  Appendices  A  and  B. 

2.4  PARAMETRIC  MODELING. 

The  basic  model  assumes  a  constant  ratio  of  specific  heats  7;  parametric  calculations  with 
different  values  for  7  can  be  conducted  using  a  modified  version  of  the  code,  where  the 
input  for  7  is  expressed  as  a  vector  rather  than  as  a  scalar,  to  study  the  sensitivity  of 
venting  to  changes  in  the  specific  heat  ratio.  Similarly,  other  parametric  calculations  with 
variations  in  fluidized  mass  fraction  /,  wall  impedance  ‘w’,  heat  diffusivity  ‘a’,  convection 
coefficient  ‘u’,  and  yield  ‘W’  may  be  performed  with  appropriately  modified  versions  of  the 
code  to  obtain  sensitivities,  to  variations  in  each  quantity,  of  the  total  vented  mass,  the 
rate  at  which  such  venting  might  occur  and  the  duration  of  such  venting. 
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SECTION  3 

THE  MATVENT  CODE 


The  computer  code  embodying  the  venting  model  for  responding  or  non-responding  cham¬ 
bers  has  been  named  MATVENT  since  the  code  was  written  in  the  high-level  language 
called  MATLAB  4.0  (Reference  3).  MATLAB  (“MATrix  LABoratory”)  is  an  interactive 
system  whose  basic  data  element  is  a  matrix.  MATLAB  provides  an  easy-to-use  computing 
environment  for  high-performance  numeric  computation  and  visualization  where  problems 
and  solutions  can  be  expressed  just  as  they  are  written  mathematically,  without  traditional 
programming.  Particularly  noteworthy  with  MATLAB  is  the  ease  and  rapidity  of  generat¬ 
ing  a  variety  of  two  and  three  dimensional  graphics.  MATLAB  allows  the  importing  and 
exporting  of  programs  written  in  other  languages  such  as  FORTRAN  and  C. 

A  listing  of  the  basic  version  of  the  MATVENT  code,  embodied  in  a  script  file  (called 
an  M-  file  in  MATLAB),  named  matvent.m  is  provided  in  Appendix  B.  Several  vereions 
of  the  code  were  generated  and  used  in  the  calculations  reported  here,  one  version  for 
each  parametric  case.  No  attempt  at  optimizing  the  code(s)  have  been  made.  Further 
optimization  may  result  in  a  more  compact  code 

MATVENT  is  an  interactive  code  for  use  on  a  PC  platform,  preferably  a  486-PC  or  bet- 
The  user  is  prompted  for  inputs  pertaining  to  the  chamber  dimensions,  wall  material 
properties,  vent  area,  explosive  mass  and  yield,  hazardous  material  (or  simulant)  mass  and 
fluidized  fraction,  specific  heats  and  the  specific  heat  ratio  of  chamber  gases,  the  initial 
and  final  times  for  the  calculation,  the  error  tolerance  for  the  calculation  and  an  option 
to  see  the  progress  of  the  calculation  on  the  computer  screen.  Then  the  code  calculates 
and  sets  up  initializations  and  matrix  or  vector  arrays  for  solution  variables.  The  pres¬ 
sure  variable  at  the  beginning  of  a  timestep  is  updated  on  an  interim  basis  by  using  a 
Runge-Kutta-Fehlberg  integration  scheme  and  calling  a  function  M-file  called  intair.m 
which  contains  expressions  for  the  time  derivative  of  chamber  pressure.  Choking  is  al¬ 
lowed  for  when  proper.  It  may  be  mentioned  in  passing  here  that  this  scheme  lends  itself 
to  convenient  extension  in  problems  of  venting  from  multiple  connected  chambers.  The 
chamber  volume  is  updated  by  including  changes  over  a  timestep  due  to  thermal  expansion 
and  interim  pressure;  the  temperatui'e  decrement  during  the  timestep  due  to  heat  loss  by 
convection  and  conduction  is  then  calculated  and  a  new  interim  temperature  is  obtained. 
The  final  pressure,  density  and  temperature  at  the  end  of  a  timestep  are  (re) -calculated 
from  the  isentropic  relations  to  account  for  the  changes  in  volume  and  temperature  during 
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the  timestep.  From  the  final  values  of  pressure,  density,  and  temperature,  the  solution 
variable  histories  -  for  overpressure,  overpressure  impulse,  density,  temperature,  efflux  ve¬ 
locity,  efflux  Mach  number,  efflux  mass  flow  rate,  integrated  efflux  mass,  relative  volume 
change  (as  a  percentage  of  original  volume)  and  relative  heat  loss  (as  a  percentage  of  the 
original  blast  energy)  -  are  obtained  and  recorded  in  vector  arrays  (or  in  matrix  arrays  in 
multiple  case  input  versions  of  MATVENT).  In  addition,  if  the  final  time  is  sufficient  to 
allow  pressure  equilibration  between  the  chamber  and  the  ambient  atmosphere,  the  venting 
duration  and  the  total  mass  vented  during  the  episode  are  also  determined. 

Finally,  the  time  histories  of  the  solution  variables  are  plotted  on  the  screen,  annotated  in¬ 
teractively  with  legends  containing  the  important  input  information,  and  then  hard  copies 
of  the  plots  are  obtained  on  a  printer.  An  example  problem  and  its  solution  are  included 
with  the  code  listing. 

MATVENT  calculates  the  solutions  swiftly,  in  a  matter  of  tens  of  seconds  to  a  few  minutes 
on  a  486-PC.  Therefore,  it  is  suited  for  rapid  trade-off  calculations  to  identify  needed  design 
features  of  a  venting  experiment.  Further  detailed  information  on  a  selected  design  may 
be  obtained  through  hydrocode  calculations. 
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SECTION  4 

RESULTS  OF  CALCULATIONS 


4.1  SIMULATED  VENTING  EXPERIMENTS  AT  WES. 

A  set  of  seven  MATVENT  calculations  was  performed  for  the  following  cases,  simulating 
the  input  parameters  of  ^'Expulsion  Experiment  No.  1”  done  recently  at  the  Waterways 
Elxperiment  Station,  U.S.  Army  Corps  of  Engineers  (Reference  4).  The  rectangular  solid 
chamber  in  this  experiment  had  the  internal  dimensions  174.00cms  x  130.81cms  x  49.53cms, 
was  made  of  thick  (^  5.08  cms)  steel  plate,  and  had  a  volume  of  1.127  cubic  meters.  The 
vent  area  was  31.67  sq.cms.  The  blast  was  produced  by  a  C-4  charge  of  mass  431  gms  (0.95 
lbs.)  Glass  beads  and  other  loose  material  with  a  total  mass  of  16,116  gms  were  placed  on 
the  floor  of  the  chamber,  pai't  of  which  were  fluidized  in  the  explosion. 

4.2  MATERIAL  PROPERTIES  ASSUMPTIONS. 

It  was  assumed  in  the  calculations  that  the  C-4  explosive  released  1400  calories  per  gram 
in  the  explosion  (Reference  5).  The  chamber  walls,  made  of  steel,  were  assumed  to  have 
the  following  properties  (Reference  6): 


•  density  -  7.86  gms/cc, 

•  acoustic  speed  (compressional  waves)  -  5,943  m/sec, 

•  heat  conductivity  -  1  J/cm  •  sec  ®K, 

•  heat  capacity  -  0.449  J/gm  •  °K, 

•  coefficient  of  linear  expansion  -  0.00001206/  °K,  and 

•  thermal  diffusivity  -  0.2834  cm^/sec. 

For  air,  the  reference  properties  assumed  were: 

•  density  -  0.001225  gm/cc  at  298.15  ®K, 

•  ambient  pressure  -  1  atmosphere  (14.7  psi),  and 

•  specific  heat  at  constant  pressure  -  1.03  J/gm  °K. 


Heat  transfer  between  the  internal  gas  and  the  walls  was  governed  by: 
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•  convective  heat  transfer  coefficient  -  1.4192  J/sec  •  cm^  °K 
4.3  CALCULATIONAL  RESULTS. 

Results  of  the  set  of  calculations  are  shown  in  Figures  4-1  through  Figure  4-12.  The 
calculational  results  show  histories  of  overpressure  ratio  (i.e.,  chamber  pressure  in  excess 
of  ambient,  divided  by  the  ambient  pressure),  overpressure  impulse  ratio  (overpressure  ratio 
integrated  over  time) ,  density  ratio  (density  divided  by  ambient  density) ,  temperature  ratio 
(temperature  divided  by  ambient  temperature) ,  efflux  velocity,  efflux  Mach  number,  mass 
flow  rate,  total  mass  of  efflux,  and  where  appropriate,  relative  volume  change  (i.e.  volume 
change  at  any  time  divided  by  original  chamber  volume)  and/or  relative  heat  loss  (i.e.  heat 
energy  lost  to  the  walls  divided  by  the  explosion  energy).  The  duration  of  venting  and  the 
total  mass  vented  are  also  shown  for  variations  of  the  particular  parameter  of  each  set. 

4.3.1  Effects  of  Varying  7;  Non-Responding  Chamber. 

Figures  4-1  and  4-2  show  the  effects  of  varying  7  for  the  case  of  a  rigid,  insulated  chamber 
with  the  assumption  that  two-thirds  of  the  mass  of  beads  was  fluidized.  The  available  ex¬ 
perimental  data  from  “Expulsion  Experiment  No.  1”  are  also  plotted  for  comparison.  The 
overpressure  ratio  (Figure  4-la)  and  overpressure  impulse  ratio  (Figure  4-lb)  comparisons 
with  the  experiment  show  very  poor  agreement  over  the  entire  range  of  7.  A  value  of  1.35 
for  7  produces  an  initial  overpressure  closest  to  the  experimental  overpressure.  Therefore, 
7  =  1.35  is  the  adopted  standard  for  the  rest  of  the  parametric  study. 

Figure  4-lc  demonstrates  that  the  common  initial  value  of  the  chamber  gas  density  de¬ 
creases  at  progressively  slower  rates  in  time  as  the  assumed  value  for  7  is  decreased,  so 
that  at  the  end  of  a  venting  episode,  lower  assumed  values  for  7  result  in  higher  residual 
densities  in  the  chamber.  The  chamber  temperatures  in  Figure  4-ld,  show  a  similar  be¬ 
havior  except  that  the  initial  temperatures  are  higher  with  the  larger  values  for  7.  The 
values  for  efflux  velocities  at  any  time.  Figure  4-le  are  generally  greater  for  larger  7  but 
taper  off  rapidly  towards  the  end  of  venting.  Figure  4-lf  shows  that  the  duration  of  choked 
flow  from  the  vent  is  reduced  for  7  <  1.2.  Figure  4-lg  indicates  that  the  efflux  mass  flow 
rates  up  to  approximately  half  of  the  total  vent  time  increase  with  7,  but  at  later  times, 
approach  common  values  regardless  of  7.  The  total  vented  mass  up  to  any  time,  shown 
in  Figure  4-lh  increases  with  7.  Figure  4-lh  shows  that  the  total  efflux  mass  decreases 
(moderately)  as  7  decreases  from  1.4  to  1.1.  Figure  4-2a  shows  that  the  venting  duration 
is  not  a  strong  function  of  7.  However,  Figure  4-2b  indicates  that  the  total  vented  mass 
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increases  moderately  as  7  is  increased. 

4.3.2  Effects  of  Varying  the  Fluidized  Mass  Fraction;  Non-Responding  Chamber. 

Figures  4-3  and  4-4  show  the  influence  of  varying  the  fluidized  mass  fraction,  /,  in  a  non¬ 
responding  chamber.  The  mass  fractions  used  were  2/3,  1/3,  1/10,  1/20  and  0;  thus,  there 
was  a  progression  from  two-thirds  of  the  loose  material  being  fluidized  to  only  air  plus  the 
mass  of  the  explosive.  As  /  is  increased,  increases  are  shown  at  any  time  in  the  chamber 
overpressure  (Figure  4-3a),  the  overpressure  impulse  (Figure  4-3b),  the  density  (Figure  4- 
3c),  the  duration  of  choked  flow  (Figure  4-3f),  the  efflux  mass  flow  rate  (Figure  4-3g), 
the  total  efflux  mass  (Figure  4-3h),  the  total  duration  of  venting  (Figure  4-4a),  and  the 
total  amount  of  vented  mass  (Figure  4-4b).  The  temperature  decreases  with  increasing 
/  (Figure  4-3d).  The  efflux  velocity  (Figure  4-3e)  shows  lower  but  longer  lasting  values 
with  increases  in  /.  Overpressure  and  overpressure-impulse  comparisons  with  data  from 
“Expulsion  Experiment  No.  1”  are  still  poor  for  any  mass  fraction  /. 

4.3.3  Effects  of  Varying  the  Wall  Impedance;  Thermally  Non-Responding  Chamber. 

Figures  4-5  and  4-6  show  the  effects  of  varying  the  characteristic  acoustic  impedance, 
‘w’  of  the  wall  material  for  a  thermally  non-responding  chamber  with  the  fluidized  mass 
fraction  assumed  as  f  —  2/3.  The  impedance  was  varied  from  that  of  steel,  w,  to  O.lw. 
Starting  with  the  same  initial  values  for  each  quantity,  the  decay  rates  of  overpressure 
(Figure  4-5a),  density  (Figure  4-5c),  temperature  (Figure  4-5d)  and  mass  flow  rate  (Fig¬ 
ure  4-5g)  increase  moderately  as  the  impedance  is  decreased  substantially.  Figure  4-5b 
indicates  that  the  overpressure  impulse  is  higher  with  higher  wall  impedance.  Figure  4-5e 
demonstrates  that  the  efflux  velocity  history  is  not  strongly  dependent  on  wall  impedance. 
Figure  4-5f  shows  that  the  duration  of  choked  flow  decreases  moderately  with  decrease  in 
impedance.  Figure  4-5h  shows  that  the  total  efflux  mass  decreases  with  decreasing  wall 
impedance  although  the  variation  is  not  large.  The  venting  duration  decreases  with  in¬ 
creasing  impedance  (Figure  4-6a),  whereas  the  total  vented  mass  increases  (Figure  4-6b). 
The  relative  volume  change  of  the  chamber  (defined  as  change  in  volume  divided  by  the 
original  volume),  Figure  4-5i,  can  be  quite  large  as  the  impedance  is  decreased.  While 
Fig  4-5i  shows  the  relative  volume  change  going  from  14  to  96  percent  as  the  impedance 
is  decreased  by  a  factor  of  ten,  it  has  to  be  borne  in  mind  that  non-linear  effects,  not 
included  in  the  current  analysis,  will  moderate  very  large  deformations  to  lower  values. 
Comparisons  with  experimental  overpressure  and  overpressure  impulse  data  are  still  poor 
(Figures  4-5a  and  4-5b). 
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4.3.4  Thermal  Response  By  Conduction  Only;  Mechanically  Non-Responding  Chamber. 

Figure  4-7  demonstrates  the  effects  of  assuming  that  the  chamber,  although  mechanically 
non-responding,  responds  by  thermal  conduction  and  also  expansion  due  to  heating.  The 
fluidization  fraction  assumed  here  is  /  =  2/3.  The  thermal  diffusivities  employed  were 
1,  1/4,  1/16,  1/100  and  0  times  the  value  for  steel,  a.  The  results  show  a  large  effect 
of  thermal  response  on  the  overpressure  (Figure  4-7a),  overpressure  impulse  (Figure  4- 
7b),  density  (Figure  4-7c),  temperature  (Figure  4-7d),  efflux  velocity  (Figure  4-7e),  mass 
flow  rate  (Figure  4-7g),  and  total  efflux  mass  (Figure  4-7h),  which  decrease  greatly  as  the 
diffusivity  is  increased  from  0  to  a.  Choking  durations  (Figure  4-7f),  as  well  as  venting 
durations  and  total  vented  mass  are  also  greatly  affected  by  increases  in  diffusivity  from 
0  to  a.  The  overpressure  and  overpressure  impulse  data  from  “Expulsion  Experiment 
No.  1”  fall  within  a  band  between  0.01a  and  0.0625a,  and  follow  the  trends  reasonably 
well;  however,  since  that  range  of  diffusivities  is  unrealistic  for  steel,  it  is  believed  that 
conduction  alone  cannot  account  for  energy  losses  due  to  heat  transfer.  A  slight  volume 
decreeise  of  the  chamber  due  to  thermal  expansion  and  a  large  amount  of  heat  loss  to  the 
walls  are  indicated  by  Figures.  4-7i  and  4-7j,  respectively. 

4.3.5  Thermal  Response  By  Convection  and  Conduction;  Mechanically  Non-Responding 
Chamber. 

Results  shown  in  Figure  4-8  assume  that  the  heat  loss  is  due  to  both  conduction  and 
convection.  The  chamber  is  assumed  to  be  mechanically  non-responding.  The  coefficients 
of  convection  employed  were  0,  u/10,  u/4,  u/2  and  u,  where  u  =  1.42  J/sec  •  cm^  •  °K; 
the  fluidization  fraction  assumed  is  /  =  2/3.  The  plots  for  time  histories  of  overpressure 
(Figure  4-8a),  overpressure  impulse  (Figure  4-8b),  density  (Figure  4-8c),  temperature  (Fig¬ 
ure  4-8d),  efflux  velocity  (Figure  4-8e),  mass  flow  rate  (Figure  4-8g),  and  total  efflux  mass 
(Figure  4-8h)  show  substantial  decreases  in  those  quantities  when  the  assumed  convection 
coefficient  is  increased  in  value.  Choking  flow  durations  (Figure  4-8f)  also  show  a  decrease 
with  increase  of  convection  coefficient.  The  relative  volume  change,  due  to  thermal  ex¬ 
pansion  only  (Figure  4-8i)  appears  to  be  of  minor  significance,  whereas  the  increase  in 
relative  heat  loss  (Figure  4-8j)  with  increasing  convection  coefficient  can  be  quite  large 
(>  40%  for  the  nominal  case  with  u  =  1.42  J/sec  •  cm^  •  "K).  In  particular,  the  plots  show 
again  that  the  total  efflux  mass  is  a  strong  function  of  thermal  response,  with  the  efflux 
mciss  decreasing  significantly  as  the  convection  coefficient  is  increased.  For  the  nominal 
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value  of  1.42  J/sec  •  cm^  °K,  the  overpressure  and  overpressure  impulse  comparisons  with 
experiment,  Figures  4-8a  and  4-8b,  are  quite  good.  The  total  efflux  mass  for  this  case, 
Figure  4-8h,  is  889  gms,  and  the  bead  mass  vented  will  be  761  gms  (i.e.,  /  x  889  gms), 
which  may  be  compared  to  a  mass  loss  of  1081  gms  (which  might  be  an  upper  bound  since 
it  included  vaporized  material  cind  escaped  dust  of  unknown  amount,  Reference  7)  in  the 
experiment. 

4.3.6  Fully  Responding  Chamber. 

Figures  4-9  and  4-10  show  the  results  for  a  chamber  which  responds  both  mechanically 
and  thermally.  It  is  assumed  here  also  that  the  fluidization  fraction  /  =  2/3.  The  con¬ 
vection  coefficient  was  varied  from  u  to  zero  as  in  case  4.3.5.  The  results  show  the  same 
trends  as  case  4.3.5.  With  the  convection  coefficient  assumed  to  have  the  nominal  value  u, 
the  overpressure  ratio  (Figure  4-9a)  and  overpressure  impulse  ratio  (Figure  4-9b)  compare 
well  with  experiment.  The  final  value  of  the  overpressure  impulse  (0.65  psi-sec)  is  slightly 
lower  than  that  for  case  4.3.5  (0.678  psi-sec);  similarly  the  total  efflux  mass  (865  gms)  is 
also  somewhat  lower  than  in  case  4.3.5  (889  gms).  The  total  efflux  mass  for  the  fluidized 
beads  is  741  gms  on  a  proportional  basis  (i.e.,  /  x  865  gms)  which  is  consistent  with  the 
experimental  data.  The  total  venting  time  (Figure  4-lOa)  and  the  total  vented  mass  (Fig¬ 
ure  4-lOb)  decrease,  as  in  case  4.3.5,  as  the  convection  coefficient  is  increased.  The  relative 
volume  change  (which  includes  both  mechanical  and  thermal  components).  Figure  4-9i, 
is  seen  to  be  in  the  range  between  one  percent  for  the  highest  value  u  of  the  convection 
coefficient  to  14%  for  u  =  0.  The  relative  heat  loss,  Figure  4-9j,  is  seen  to  range  from  53% 
for  the  case  with  convection  coefficient  u  to  0  for  a  thermally  non-responding  case  (u  =  0). 

4.3.7  Scaling  Issues. 

Figures  4-11  and  4-12  depict  the  results  obtained  when  different  yields  and  corresponding 
“yield-scaled”  chambers  and  chamber  contents  are  assumed.  These  calculations  employed 
a  fluidized  fraction  /  =  2/3,  thermal  diffusivity  a  =  0.2834  cm^/sec,  convection  coefficient 
u  =  1.42  J/sec  •  cm^  •  °K,  and  wall  impedance  w  =  4-67- 10®  gm/cm^-  sec.  The  parametric 
yields  used  were  multiples  of  1,  8,  64,  343,  and  1,000  times  the  431  gms  of  C-4  used  in  the 
WES  “Expulsion  Experiment  No.  1”.  Corresponding  chamber  volumes,  explosive  masses 
and  hazardous  material  masses  also  had  the  same  multiples.  Vent  areas  were  2/3  power 
multiples  while  linear  dimensions  were  1/3  power  multiples.  The  solution  variables  were 
also  appropriately  scaled.  Time  in  each  case  was  scaled  by  the  inverse  1/3  power  of  the 
yield  multiple.  If  yield  scaling  is  valid,  a  single  curve  should  represent  the  time  history 
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of  any  particular  solution  variable  on  a  scaled  basis  for  all  yields.  This  does  not  appear 
to  be  the  case  since  in  a  responding  chamber,  the  losses  of  energy  due  to  mechanical  and 
thermal  processes  do  not  scale  with  yield.  In  general,  the  yield-scaled  venting  duration 
and  total  mass  vented  increase  eis  the  yield  increases,  as  shown  by  Figures  4-12a  and  4-12b, 
respectively. 
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Overpressure  Impulse  Ratio  Overpressure  Ratio 


Figure  4- 1 .  Blast  venting  —  effect  of  varying  gamma. 
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Figure  4-1.  Blast  venting  —  effect  of  varying  gamma  (Continued). 
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Efflux  Velocity,  cm/sec  (10^) 


Figure  4-1.  Blast  venting  —  effect  of  vaiying  gamma  (Continued). 
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Total  Vented  Mass,  gms  (io^)  ^  Total  Venting  Time,  sec 


1,25 

gamma 


a.)  Total  venting  time. 


1.25 

gamma 


b.)  Total  vented  mass. 


Figure  4-2.  Dependence  of  vent  duration  and  vented  mass  on  gamma. 
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b.)  Overpressure  impulse  ratio  history. 


Figure  4-3.  Blast  venting  —  effect. of  mass  loading. 
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Temperature  Ratio  Density  Ratio 


Figure  4-3.  Blast  venting  —  effect  of  mass  loading  (Continued). 
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Mach  Number  of  Efflux 


Time,  sec 


f.)  Mach  number  history. 

Figure  4-3.  Blast  venting  —  effect  of  mass  loading  (Continued). 
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Figure  4-3.  Blast  venting  —  effect  of  mass  loading  (Continued). 
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bead  fraction  fluidized 


b.)  Total  vented  mass. 


Figure  4-4.  Dependence  of  vent  duration  and  vented  mass  on  fluidized  bead  fraction. 


a.)  Overpressure  ratio  history. 


b.)  Overpressure  impulse  ratio  history. 

Figure  4-5.  Blast  venting  —  effect  of  wall  compliance. 
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Figure  4-5.  Blast  venting  —  effect  of  wall  compliance  (Continued). 
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Figure  4-5.  Blast  venting  —  effect  of  wall  compliance  (Continued). 
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g.)  Mass  flow  rate  history. 


Figure  4-5.  Blast  venting  —  effect  of  wall  compliance  (Continued). 
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Figure  4-5.  Blast  venting  —  effect  of  wall  compliance  (Continued). 
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normalized  wall  impedance 


b.)  Total  vented  mass. 


Figure  4-6.  Dependence  of  vent  duration  and  vented  mass  on 


Figure  4-7.  Blast  venting  —  heat  loss  by  conduction  only. 
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c.)  Density  ratio  history. 
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d.)  Temperature  ratio  history. 


Figure  4-7.  Blast  venting  —  heat  loss  by  conduction  only  (Continued). 
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Figure  4-7.  Blast  venting  —  heat  loss  by  conduction  only  (Continued). 
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Figure  4-7.  Blast  venting  —  heat  loss  by  conduction  only  (Continued). 


Figure  4-7.  Blast  venting  —  heat  loss  by  conduction  only  (Continued). 
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b.)  Overpressure  impulse  ratio  history. 


Figure  4-8.  Blast  venting  —  heat  loss  by  convection  and  conduction. 
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Figure  4-8.  Blast  venting  —  heat  loss  by  convection  and  conduction  (Continued). 
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Figure  4-8. 


Blast  venting  —  heat  loss  by  convection  and 


conduction  (Continued) 
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Figure  4-8.  Blast  venting  —  heat  loss  by  convection  and  conduction  (Continued). 
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Figure  4-8.  Blast  venting  —  heat  loss  by  convection  and  conduction  (Continued). 
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c.)  Density  ratio  history. 


d.)  Temperature  ratio  history. 


Figure  4-9.  Blast  venting  —  fully  responding  (Continued). 
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Mach  Number  of  Efflux  Efflux  Velocity,  cmfaec  (10^) 


f.)  Mach  number  history. 


Figure  4-9.  Blast  venting  —  fully  responding  (Continued). 
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Figure  4-9.  Blast  venting  —  fully  responding  (Continued). 
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i.)  Relative  volume  change  history. 


Figure  4-9.  Blast  venting  —  fully  responding  (Continued). 
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b.)  Total  vented  mass. 


Figure  4-10.  Dependence  of  vent  duration  and  vented  mass  on 
convection  coefficient. 
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Scaled  Overpressure  Impulse  Scaled  Overpressure  Ratio 


a.)  Overpressure  ratio  history. 
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b.)  Overpressure  impulse  ratio. 


Figure  4-11.  Scaled  yield  effects  on  venting  —  fully  responding. 
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Scaled  Temperature  Ratio  Scaled  Density  Ratio 


c.)  Density  ratio  history. 


Scaled  Time 


d.)  Temperature  ratio  history. 


Figure  4-11.  Scaled  yield  effects  on  venting  —  fully  responding  (Continued). 
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Figure  4-11.  Scaled  yield  effects  on  venting  —  fully  responding  (Continued). 
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g.)  Mass  flow  rate  history. 


h.)  Total  efflux  mass  history. 


Figure  4-11.  Scaled  yield  effects  on  venting  —  fully  responding  (Continued) 


i.)  Relative  volume  change  history. 
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b.)  Vented  mass. 
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SECTION  5 
DISCUSSION 


The  primary  results  obtained  in  the  parametric  studies  of  Section  4  are  summarized  in 
Table  5-1.  Venting  is  most  significantly  affected  by  two  factors,  namely  the  fluidized 
fraction  and  the  convection  coefficient.  Variations  in  7  and  the  wall  impedance  appear  to 
be  less  important  than  the  properties  related  to  thermal  effects,  such  as  the  convection 
coefficient  and  the  diffusivity,  in  their  effect  on  venting  durations  and  vented  masses. 
Conversion  of  the  results  of  small  scale  experiments  to  large  scale  by  simple  yield-scaling 
becomes  more  complicated  because  of  the  non-scaling  energy  losses  due  to  thermal  effects. 
In  particular,  any  predictions  for  venting  durations  and  total  vented  masses  for  a  large 
event  based  on  yield-scaling  of  a  small  event  are  likely  to  be  underpredictions,  as  shown 
by  Figure  4-12. 

The  convection  coefficient  u  is  a  quantity  which  depends  in  general  on  the  geometry  of  the 
chambers,  the  physical  and  thermal  properties  of  both  chamber  gases  and  chamber  walls, 
the  conditions  of  the  wall  surfaces,  and  the  fluid  mechanical  regime  in  the  chamber.  It  is  not 
a  well-established  quantity  for  chamber  explosions,  although  considerable  effort  has  gone 
into  its  determination  in  other  practical  engineering  applications  (Reference  8).  In  forced 
convection  problems,  u  is  related  to  the  Nusselt  number.  The  range  of  values  for  convection 
coefficients  is  also  quite  large,  going  from  0.0003  J/sec  •  cm^  •  °K  for  heating  (or  cooling) 
of  air  to  28.4  J/sec  •  cm^  •  °K  for  dropwise  condensation  of  steam  (Reference  9).  A  value  of 
9.23  J/sec  •  cm^  •  °K  has  been  suggested  for  the  flow  of  liquid  metals  in  pipes  (Reference  2). 
The  value  1.42  J/sec  ■  cm^  •  °K  appears  to  give  MATVENT  results  in  agreement  with  data 
from  “Explusion  Experiment  No.  1”  for  a  presumed  fluidization  fraction  of  2/3.  At  present 
no  experimental  data  are  available  to  substantiate  independently  the  selection  of  such  a 
value.  Other  values  (lower)  of  fluidization  fraction  and  convection  coefficient  may  also  give 
results  in  agreement  with  experimental  data.  The  convection  coefficient  depends  on  the 
mass  flow  rate  in  forced  flow  pi'oblems.  Therefore,  it  is  very  likely  that  the  convection 
coefficient  for  blast  in  responding  chambers  is  dependent  on  the  density  and  fluidized 
fraction  of  the  hazardous  material.  For  very  small  amounts  of  fluidized  material,  when 
the  density  of  chamber  gases  approach  the  density  of  air,  it  is  likely  that  the  convection 
coefficient  is  two  orders  of  magnitude  smaller.  There  is  a  need  for  searching  the  literature 
and  developing  a  formulation  for  the  convection  coefficient  for  problems  of  explosions  in 
venting  chambers. 


55 


The  loose  solid  material  in  the  chamber  has  been  addressed  in  this  report  within  the  sim¬ 
ple  “dense  -gas”  approximation  framework,  wherein  the  solid  material,  explosion  gases  and 
chamber  air  behave  as  if  the  mixture  constituted  a  single  dense  gas,  with  no  drag  or  heat 
transfer  between  the  particle  and  gas  phases.  The  next  level  of  analytical  sophistication 
would  be  to  employ  the  theory  of  interacting  continua,  where  the  solid  and  gaseous  phases 
are  treated  as  co-existing  continua  occupying  the  same  volume,  but  capable  of  exchanging 
momentum  by  means  of  drag  forces  and  energy  by  means  of  heat  transfer.  In  reality, 
however,  the  particles  are  discrete,  and  particle-particle  collision  processes  also  will  have 
to  be  taken  into  account,  especially  when  the  mass  of  the  solid  particles  becomes  large,  in 
the  exchange  of  momentum  and  energy.  If  inelastic  particle  collisions  are  included,  energy 
losses  due  to  plastic  deformations  have  to  be  taken  into  account.  Further  complications  are 
presented  by  the  possibility  of  particle  agglomerations,  breakups,  and  subsequent  size  re¬ 
distributions,  all  within  a  highly  turbulent  confined  flow  produced  by  an  explosion.  Phase 
changes  from  solid  to  liquid  to  gas  can  occur  in  regions  where  the  temperatures  are  high 
enough  and  durations  long  enough,  while  condensation  and  solidification  might  occur  in 
cooler  regions.  Estimations  of  heat  transfer  even  to  a  single  spherical  particle  of  known 
thermal  conductivity,  density  and  specific  heat  capacity  can  be  difficult  because  of  uncer¬ 
tainties  associated  with  surface  characteristics  (represented  in  the  convection  coefficient). 
For  example,  in  the  absence  of  any  surface  resistace  to  heat  transfer,  a  glass  sphere  of 
0.01  cm  diameter  (k  =0.85  J/sec-  cm  •  °K  and  a  =  0.0046  cm^/sec,  approximately,  for 
common  glass)  attains  temperature  equilibrium  with  a  surrounding  hot  air  medium  in  ap¬ 
proximately  5  msec;  however  the  same  sphere  will  not  heat  up  at  all  if  the  surface  were  to 
be  insulated.  These  uncertainties  are  not  addressed  by  the  simple  model  presented. 

The  MATVENT  model  has  also  not  considered  the  following  factors:  (a)  loss  of  energy  due 
to  heating  and  vaporization  of  the  hazardous  material,  (b)  heat  addition  to  the  chamber 
gases  due  to  any  exothermic  reactions  in  the  material  being  vaporized,  (c)  heat  additions 
due  to  condensation  of  hazardous  material  vapors  near  the  (cooler)  wall  surfaces,  (d)  the 
effect  of  increzised  wall  ai'ea  accompying  changes  in  the  chamber  volume,  (e)  the  effect  of 
any  changes  in  the  vent  area  during  venting,  and  (f)  the  effect  produced  by  constriction  of 
the  jet  of  venting  gases  at  the  orifice. 

For  high-yield  explosions  (a),  (b)  and  (c)  are  probably  minor  effects.  In  the  case  of  (d),  a 
separate  calculation  (not  reported  here)  for  fully  responding  chambers  in  which  the  wall 
area  was  allowed  to  change  proportionally  as  the  two-thirds  power  of  the  chamber  volume, 
failed  to  show  any  significant  differences  from  the  “fully  responding”  case  (Section  4.3.6). 
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Table  5-1.  Summary  of  parameter  variation  effects  on  venting. 


Quantity 

increasing 

11^1112^^^^1111 

increasing 

w 

increasing 

a 

increasing  1 
u 

Overpressure 

increases 

increases 

initially 
same, then 
increases 

decreases 

decreases 

Overpr  Impulse 

increases 

increases 

increases 

decreases 

decreases 

Density 

decreases 

increases 

initially 

same.then 

increases 

decreases 

decreases 

Temperature 

increase 

then 

decrease 

decreases 

initially 

same,then 

increases 

decreases 

decreases 

Efflux  Velocity 

increases 

increase 

then 

decrease 

initially 

same.then 

increase/ 

decrease 

decreases 

decreases 

Choking 

Duration 

increases 

increases 

increases 

decreases 

decreases 

Mass  flow  rate 

increases 

increases 

increases 

decreases 

decreases 

Total  Efflux  mass 

increases 

increases 

increases 

decreases 

decreases 

Volume  change 

— 

— 

decreases 

decreases 

decreases 

Heat  loss 

- 

— 

— 

increases 

increases 

Venting  duration 

little  effect 

increases 

decreases 

decreases 

decreases  | 

Total  Vented 

mass 

increases 

increases 

increases 

decreases 

decreases  | 

Vent  area  changes  during  efflux  may  be  significant  in  estimating  vented  masses  and  venting 
durations;  a  reasonable  model  for  such  changes  needs  to  be  developed.  Jet  constriction 
effects,  which  depend  upon  vent  pipe  length  and  geometry,  can  be  taken  into  account  by 
means  of  an  “orifice  coefficient”  such  as  that  for  a  Borda  tube.  It  is  possible  to  incorporate 
all  of  the  foregoing  effects  into  future  versions  of  MATVENT. 
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SECTION  6 

CONCLUDING  REMARKS 


A  simple  PC-based  code  named  MATVENT  has  been  developed,  described  and  employed 
to  carry  out  parametric  calculations  on  the  effects  of  varying  7,  varying  the  fluidization 
fraction,  varying  the  chamber  wall  impedance,  varying  the  chamber  wall  diffisivity  and 
varying  the  convection  coefficient  in  several  separate  studies.  The  response  models  provide 
an  approximate,  but  convenient,  method  in  assessing  the  relative  importance  of  factors 
which  affect  venting.  The  approximate  nature  of  the  models  is  consistent  with  the  current 
state  of  knowledge  on  venting  from  bleist  in  chambers.  The  effects  of  chamber  response  on 
yield-  scaling  of  venting  durations  and  vented  masses  is  also  examined  with  MATVENT 
runs.  It  is  shown  that  yield-scaling  based  on  small  scale  experiments  may  lead  to  signifi¬ 
cant  underpredictions  of  venting  durations  and  vented  masses  in  a  large  scale  event.  The 
parametric  calculations  show  that  attention  must  be  paid  to  thermal  response,  and  to  a 
lesser  extent,  the  mechanical  response,  of  the  chamber  wall  materials  in  designing  small 
scale  experiments  on  internal  blast  and  venting.  The  chamber  materials  must  be  carefully 
chosen  and  the  internal  surfaces  of  the  chamber  must  be  carefully  designed  to  achieve  the 
purpose  of  the  experiments. 

MATVENT  provides  a  tool  for  quick  (a  minute  or  so  per  calculation  on  a  486  PC  platform) 
trade-off  studies  in  experiment  design.  MATVENT  can  also  be  useful  when  conventional 
scaling  rules  become  uncertain  in  the  calculation  of  potential  full  scale  scenarios. 

Further,  MATVENT  can  be  extended  to  treat  cases  of  venting  from  multiple  intercon¬ 
nected  chambers  in  one  of  which  a  blast  might  be  set  off.  Further  research  needs  exist 
in  establishing  a  formulation  for  convection  coefficients  under  chamber  blast  conditions, 
and  in  particular,  the  dependence  of  convection  coefficients  on  the  properties  of  chamber 
gases.  Extensions  to  the  current  MATVENT  model  will  be  helpful  in  assessing  the  relative 
sensitivity  of  venting  durations  and  vented  masses  to  factors  such  as  vent  area  changes 
during  venting,  jet  constriction  effects  and  vaporization,  condensation  as  well  as  reaction 
characteristics  of  stored  hazardous  materials. 
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APPENDIX  A 

MODELING  EQUATIONS 


W  =  Pw  •  Cw 
^  =  u/{py,  ■  C  •  ® 

POc  =  {pOaVd  +  Afci}/Vb 
■Ro  =  POp/POa 
COap  =  {7  -POo/POpj^'^ 

^Op=[l+{(7-l)-W/(yo-pOa)}] 

Po  =  POa  •  a^Op 

a;  - 

•‘-cp  —  •‘'Op 

COp  —  COap  "  yJ^Og 

CQcr  =  COg  •  V^2/(7-f  1) 

eo=[i  +  'r-w/{Cp-pog-Vo^Toa)] 

To  =  Toa  •  ©0 

©cp  ^  ^cp  *  ©0 
Xg  =  PgIpOa 

9,  =  T,(()/T„ 

A)  =  { (C/C,)  ■  (p„/pa,)  ■  (ARo/Vo)  ■  ^a} 

A  =  ARo‘Poa/w 

B  =  ARq  •  Toa  •  €  •  y/a 

X={{7  +  l)/2r^‘-'-‘' 

A)ar  =  POg  '  [2/(7  + 

Uo  —  Cocr 

mo  =  ao  ’  pocr  ■  Uq 

Mq  =  1 
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if  iCOfl  <  X, 


rho—aa-  \/(7+  l)/(7-  1)  •  POg  •  cocr  • 

Uo  =  7ho/ (ao  •  POer) 

Mo  =  N/2/h  -  1)  ■  -  l}*'" 


-2/-r  -(7+l)/7-|  1/2 

Os  "  ^Og  J 


Fluid  Miechanical  Model-Quasi  steady  trcament  for  the  rate  of  change  of  pres¬ 
sure  in  the  chamber: 

Xg{t)  =  X(t) 

if  ©s(^)  >  ^cgy 

%(t)  =  -7{(7  +  .  c,(t)  ■  (A/V)  ■ 

if  @,(t)  <  e^, 

Xp  =  0 

Where,  if  x{t)  >  x, 

A{t)  —  ao 

if  ^(t)  <  X, 

A{t)  =  ao  •  V2/(7  -  1)  •  {(7  +  1)/2}<^+i)/2(t'-i)  •  -  1 

Cp  =  cop  •  (xp/rrop)^^"^^^'^ 

The  mechanical  response  model  consists  of  two  parts: 

1.  Reduction  of  chamber  volume  due  to  thermal  expansion  of  the  “penetration  thickness” 
ya  •  t  during  a  time  interval  St: 


0(£)  = 

SVte  =  -B  •  \/i  •  (0pi  -  l)  •  (1  -  0)  •  5£ 

2.  Change  in  chamber  volume  during  the  interval  St  due  to  pressure  change: 

SVpr{t)  =  A  ■  (xg  —  l)  ■  St 
The  net  volume  change  over  interval  St  is: 

SV{t)=SVte{t)  +  SVprit) 
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The  chamber  volume  at  time  t  is  then: 

V{t)  =  V{t  -  6t)  +  6V{t) 

Heat  loss  modeling  is  as  follows: 

The  thermal  response  of  the  chamber  due  to  heat  transfer  at  the  walls  can  be  modeled  as 
being  due  (a)  either  to  conduction  only,  or  (b)  due  to  both  convection  and  conduction. 

The  temperature  ratio  change,  ^0p(t),  in  the  chamber  during  a  time  interval  6t  due  to 
heat  loss  to  the  walls  is  calculated  as  follows: 

©rW  =  {©si(<)  +  ©p(t-^i)}/2 

For  t  >  0,  if  0“"(t)  >  0cp, 

(a)  if  it  is  assumed  that  heat  loss  is  due  to  conduction  only,  then 

8Qg{t)  =  -/3o  •  .  {Qgi  -  l)/2v/i)  +  {Qgi  -  0“")  •  sTt] 

(b)  if  the  assumption  of  Newtonian  heating  of  the  chamber  walls  with  surface 
conductance  (i.e.  convection  and  conduction)  is  made,  then 

6Qg{t)  =  ARo/iPg  -V-Cp)}-  (0“"  -  1) .  O 

If  Qg^{t)  <  Qcgi  or  if  t  <  0, 

=  0 

The  net  temperature  ratio  change  A0p  due  to  heat  loss  up  to  time  t  is 

ABgit)  =  A0p(t  -  5t)  +  SSgit) 

The  net  temperature  change  as  a  fraction  of  the  initial  (maximum)  temperature  is: 

DQg{t)  =  AQg{t)/eo 
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The  chamber  and  efflux  variables,  adjusted  for  mechanical  and  thermal  response  of 
the  chamber,  can  be  calculated  as  follows: 

(a)  the  adjusted  chamber  pressure  ratio  x 

x(t)  =  I,  .  [{(K  -  6V^)lvy  ■  {max(e,  +  6%, 

(b)  the  adjusted  overpressure  ratio  Ox 

Ox(t)  =  x(t)  —  1 

(c)  the  overpressure  impulse  ratio  OI 

01  (t)  =  01{t  —  +  Ox  •  6t 

(d)  the  adjusted  chamber  density  ratio  R 

R{t)  =  Ro  •  {x/xog)^/'^ 

(e)  the  adjusted  critical  density  Pcr 

l>or{t)  =  Po,  ■  (R/Ro)  ■  {2/(7  + 1)}'/'’-*' 

(f)  the  adjusted  chamber  temperature  ratio  0 

0(t)  =  00  •  {x/xog)^'^~'^'^^'^ 


(g)  the  sound  speed  in  the  chamber  c 

—  Cpg  • 


(h)  the  critical  sound  speed  Ccr 
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the  velocity  of  efflux  U,  the  mass  flow  rate  rh,  and  the  efflux  Mach  number  M: 


U{t)=c^ 
m(t)  =  ao-  Per' U 

M{t)  =  1 


if  a:  <  X, 


m(t)  =  ao  •  v/(7  +  l)/(7  -  1)  -poa'R  cer 
U{t)  =m/{ao'  Per) 

M(t)  =  ^{2/(',-l)}-(a;-i/T-l) 


the  meiss  flow  increment  6m  and  the  total  meiss  outflow  by  time  t 

6m{t)  =  m-St 

=  rntot  {t  -  6£)  +  6m 


APPENDIX  B 


LISTING  OF  THE  CODE  MATVENT 

(With  an  Example  Run) 


%matvent.m 

% _ 


%AN  INTERNAL  AIRBLAST  CODE  FOR  VENTING  FROM 
%OR  NON-RESPONDING  SINGLE  CHAMBER. 

% 


A  RESPONDING 


^DEVELOPED  BY  —DEV  S.  SRINIVASA 

%LOGICON  RDA,  6053  W.  CENTURY  BLVD,  LOS  ANGELES,  CA  90045 
%PH:  (310)  645-1122  ext  338;  FAX(310)  645-2553 
% _ _ _ 


%Set  plotting  indices  as  elenic7its  of  a  1X5  matrix  iji  which; 

%elcment  1  is  the  index  for  overpressure  and  ovpr  impulse  ratio  plots 
%element  2  is  the  index  for  density  and  temperature  ixitio  plots 
%element  3  is  the  index  for  efflux  velocity  and  Mach  number  plots 
%element  4  is  the  index  for  mass  flow  rate  and  total  efflux  mass  plots 
%clcment  5  is  the  index  for  relative  voluttie  ehange  and  heat  loss  plots 


kplt-inpnt(’plot  indices  as  a  1x5  inatrix’);%for  example, type[l  111  l],liit  Enter 


Vo  The  input  values  for  the  example  me  shown  after  % 

MSIM=inpiit(’Mass  of  simulant, gms’);%  16116 
frac=inpnt( Traction  fluidized’);%  0.667 
MC4=input(’Mass  of  explosive,  gms’);%  0.95*453.6 
ECON=input(’Yiel(l  per  gm  of  explosive,like  6.65el0  for  C-4’)%  6.65el0 
W=MC4*ECON;  %Explosion  yield  of  0.95  lbs  C-4 

Cst=input(’Thermal  capacity  of  chamber  wall  material  J/gm./deg  K,like  0.449e7  for  iron  or  steer);%0.449e7 

kst=input(’Tlmrinal  conductivity  of  wall  material,ergs/cm.sec.K.,like  1.0e7  for  steer);%1.0e7 
l.xpst-mput(  Coefft  of  linear  expansion,cm/cm/  K,like  1.206e-5  for  1020  steer);%1.206e-5 
diost_mput(’Density  of  wall  material, gms/cc,like  7.86  for  1020  steel ’);%7.86 
/ohote:  Set  iisk  to  0  for  thermally  non-respondiny  case 

usk=input(’Convection  coefft ,ergs/sec.cnr 2  K,like  1.419222e7  for  steel’)  %1.419222e7 
ast-kst/(rhost*Cst);%Thermal  diffusivity  of  wall  material, cnr2/sec 
flm=usk/(rhost*Cst*kst)'0.5;  %film  coefficient 

Cp  =  1.03e7;  %1.03  joules/gm  K  (CRC  Handbook),Sp.  heat  at  con.pr  for  air 
c-input(  Chamber  length,width, depth  as  a  1x3  matrix,  all  in  cms’)%[173.99  130.81  49  531 

0  lc(l,l)  lc(l,2)  lc(l,3)  %Initial  chamber  volume,  cc 
lhmd=6*V0/Arst;  %length  scale  for  chamber  (6*Vol/Area)  cm 
rhoza=0.001225  %  Ambient  air  density, g/cc 

rhozc=(rhoza*V0+MC4)/V0;  %initial  density  of  chamber  contents  w/o  beads 
rhozt0=rhozc+(frac*MSIM/V0);  %Density  of  fluid  plus  fluidised  material 
area=input  (’Area  of  vent  opening,sq.cms’);%31.66922 
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p0=1.01325e6;  %Referenc.e  pressure  in  dynes/sq.cm  (=1  atm) 
gam=input(’Assumed  gamma  for  chamber  gases’);%1.35 
TEMR=:298.15;  %Reference  temperature  25  C 
betaO=(Cst/Cp)*(rhost/rhoztO)*(Arst/VO)*sqrt(ast); 

ssO=sqrt(gam*pO/rhoztO);  %Sound  speed  in  ambient  chamber  gases  at  atm  •  pr  • 
xO=l+((gam-l)*W/(VO*pO));  ^Initial  (Max)  Pressure  ratio 
xcut=(l/xO)‘((gam-l)/gam);%cutoff  value  for  pressure 
ssmaxO=ssO*sqrt(xO);  %average  initial  sound  speed  in  chamber 
sscrit0=((2/(gam+l))‘0.5)*ssmax0;  %Initial  critical  speed  of  sound 


%Note:  set  ssch  to  a  very  large  value, say  l.e25,  for  non-responding  case 

ssch=input(’Compressional  wavespeed  in  wall  material, cm /s, like  5.9436e5  for  steel’);  %5.9436e5 
chimpd=rhost*ssch;  %Chamber  matl  impedance  gm/cm**2.sec 
cndv=(Arst*pO/chimpd);  %constant  in  volume  change  adjustment  to  pressure 
DRO=(rhoztO/rhoza);  %Max  (initial)  density  ratio 

DTO=gam*W/(Cp*rhoztO*VO);  %Temp  rise, assuming  that  all  fluid  mass  in  chmbr  participates 
TEMO— (DTO+TEMR);  %Initial  temperature(deg  K)  of  chamber  after  explosion 
TRO=TEMO/TEMR;  %InitiaI  temperature  ratio 
TRcut=xcut*TRO; 

gchok=(0.5*(gam+l))  (gam/(gam-l));  %  Choking  pressure  ratio. 
rhoztOc=((2/(gam+l))*(l/(gam-l)))*rhoztO;  %Critical  density 
if  (xO  /gchok)  >  1, 

VELO=sscritO;  %Efnux  velocity  is  same  as  critical  sound  speed 
MSDOTEO=area*rhoztOc*VELO;  %Mass  rate  of  flow 
.ME0=1; 

elseif  (xO  /gchok)  <  1, 

MSDOTE0=area*(((gam+l)/(gam-l))‘0.5)+(rhozt0*(sscrit0*((x0*(-2/gam))-(x0‘(-(gam+l)/gam)))*0.5)); 

VELO=MSDOTEO/(rhoztOc*area); 

MEO=((2/(gain-l))*(xO‘((gam-l)/gain)-l))"0.5; 

else, 

VELO=sscritO;  %Efllux  velocity  at  choking  point 
MSDOTEO=area*rhoztOc*VELO; 

ME0=1; 

end; 


%Specify  initial  time, end  time,tolcrance,trace  in  IXJ,  matrix 

TMTOL=input(’initial  time, end  time, tolerance, trace  of  calcii  as  1X4  matrix’);%[0  1  l.e-4  0] 

tO=TMTOL(l,l);tfinal=TMTOL(l,2);tol=TMTOL(l,3);trace=TMTOL(l,4); 

%  Initialization 
pow  =  1/3; 
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t  =  tO; 

hmax  =  (tfinal  - 1)/100; 
limin  =  (tfinal  -  1)/20000; 
h  =  (tfinal  -  t)/100; 

xx0=x0’;DDR0=DR0’;TTR0=TR0’;VVEL0=VEL0’;MME0=ME0’; 

X  =  x0(:);  ox=x-l; 

beta=betaO(:); 

ssmax=ssinaxO(:); 

rhozt=rlioztO(:); 

V  =  V0(:); 

TR=  TR0(:); 

DR=  DR0(:); 

VEL=VEL0(:); 

ME=  ME0(:); 

MSDOTE=MSDOTEO(:); 
tout  =  t; 

xoiit  =  X.’;  len=lengtli(xout); 
oxout=ox.’; 


for  l=l:len, 

Tls0(i)=0; 

DTls0(i)=0; 

TTMS0(i)=0; 

010(i)=0; 

DV0(i)=0; 

end; 

Tls=TlsO(:);DTls=DTl.sO(:);  TOTMASE=TTMSO(:);OI=OIO(:);DV=DVO(:); 
DRout=DR.‘  ; 

TRout=TR.’  ; 

VELout=VEL.’; 

MEout=ME.’; 

MSDout=MSDOTE.’; 

TOTMoiit=TOTMASE.’; 

OIout=OI.’: 

DVout=(DV./V0)’; 

DTout=(Tls./TRO)’; 


tail  =  tol  *  niax(norm(x,  ’inf’),  1); 
if  trace 
clc,  t,  ll,  X 
end 


%The  maw  loop 

while  (t  <  tfinal)  k  (h  >=  liinin) 
if  t  +  h  >  tfinal, 
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h  =  tfinal  -  t; 
end 

TRB=TR; 


Vo  Compute  the  slopes,  using  the  Runge-Kutta  Fehlberg  integration  scheme 


si  =  feval(’intair8’,  t,  x,  V,xxO,ssmax, area, TR,TRcut, gam);  si  =  sl(:); 

s2  =  feval(’intair8’,  t+li,  x+li*sl,  V,xxO,ssmax, area, TR,TRcut, gam);  s2  =  s2(:); 

s3  =  feval(’intair8’,  t+h/2,  x+li*(sl+s2)/4,  V,xxO,ssmax, area, TR,TRcut, gam);  s3  =  s3(:); 


%Estimatc  the  error  and  the  acceptable  error 


delta  =  norm(h*(sl  -  2*s3  +  s2)/3,’inr); 
tan  =  tol*max(norm(x,’inf’),1.0); 


%  Update  the  solution  only  if  the  error  is  acceptable 


if  delta  <  =  tan 
t  =  t  +  h; 

X  =  X  +  li*(sl  +  4*s3  +  s2)/6; 


%Reduction  of  Chamber  xiolume  due  to  thermal  expansion  of  pentrn  thekness 


dvtcon=Arst*lxpst*sqrt(ast)*TEMR;  ex=exp(-2*flin*sqrt(t)); 
d  vt =-cl  vtcon  *li  *(  t  ‘0.5)*(TR- 1  )*•(  1-ex); 


%pressure  adjustment  for  volume  change  during  timestep  h 


V=V  -|-li*ciulv*(x-l)-t-dvt;  VR=(V/V0);  DV=(V-V0)./V0; 


%  Temp  change  in  chamber  during  time  increment  h  due  to  heat  loss 


TRAVG=0..'3*(TR-|-TRB); 
if  t  >  0.0, 

if  TRAVG  >  TRcut, 


%  Newtonian  heating  only  of  chamber  wall,  with  surface  conductance 
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DTls=-h*(uslc*Arst/(rliozt*V*Ci)))*(TRAVG-l)*exp(-2*flm*t*0.5); 

elseifTRAVG<TRcut, 

DTls=0; 

else, 

DTls=0; 

end; 

elseif  t  <  0.0, 

DTls=0.0; 

else, 

DTls=0.0; 

end; 

Tls=Tls+DTls; 

DTl=Tls./TR0; 


%Adjusted  pr-essui'c  ratio 

x=real(x*((((V-h*c.iKlv*(x-l))/V)‘gain)*(((niax(TR+DTls,TRcut))/TR)*(  gam/(gam-] ))))); 

%Adjustcd  density  ratio 

DR=  real(DDRO*((x/xxO)‘(l/gani))); 

%Adjustcd  temperature  ratio 

TR=  real  ( TTRO  *((x7xx0 )  ‘  ( (gam- 1 )  /gam ) ) ); 

%sound  speed  in  the.  ehamber 

ssmax=real(s.smaxO*(x*DDRO/(xxO*DR))'0.5); 

%  “critical”  speed  oj  sound; 
sscrit  =  (sqrt(2/(gam-t-l)))*ssmax; 

%Overpressure  ratio 


ox=x-l; 


%  Overpressure  impulse  ratio 
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OI=OI+ox*h; 


YoCritical  density  at  updated  time 
rlioztc=rhozt0*((2/(gam+l))‘(l/(gam-l)))*(DR/DDR0); 


%  Velocity  and  Mass  flow  rate 
if  X  >  gchok, 

VEL=sscrit;  %Efflux  velocity  is  same  as  critical  sound  speed 

MSDOTE=area*rhoztc*VEL;  %Mass  rate  of  flow 

ME=1; 

elseif  X  <  gchok, 

MSDOTE=real(area*(((gain+l)/(ga.m-l))‘0.5)*rhoza.*(DR*(sscrit*(x'(-2/gam)-x‘(-(gam+l)/gam))'0.r)))' 

VEL=MSDOTE/(rhoztc*area); 

ME=real(((2/(gam-l))*((x)‘((gani-l)/gam)-l))*0..'5); 

else, 

VEL=sscrit;  %  Efflux  velocity  at  choking  point 
MSDOTE=area*rhoztc*VEL; 

ME=1; 
end ; 


YoMass  floxL'  increment 


DMSE  =  MSDOTE*h; 


YoMass  outflow  upto  time  t 


TOTMASE=TOTMASE  +DMSE; 


Yistore  t>ent  times  and  total  77iass  vented 


if  ox  >  =  0.01 

Vtni=t;  Tmvtd=TOTMASE; 
end; 


tout  =  [tout;  t]; 
xout  =  [xout;  X.’]; 
oxout=  [oxout;ox.’]; 
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OIout=[OIout;OI.’]; 

DRout=  [DRoiit;  DR.’]; 

TRout=  [TRout;  TR.’]; 
VELout=[VELout;  VEL’j; 

MEout=  [MEout;  ME.’]; 

MSDout=[MSDout;MSDOTE.’]; 

TOTMout=[TOTMout;TOTMASE.’]; 

DVout=[DVout;DV.’]; 

DTout=[DTout;DTl.’]; 

end  %end  of  the  ’if  delta  <  tan’  loop 


if  trace 
home,  t,  h,  X 
end 


%  Update  the  step  size 


if  delta  ~  =  0.0 

h  =  min(hmax,  0.9*h*(tan/delta)‘pow); 
end 

end;  %end  of  the  main(‘\vhile’)  loop 


if  (t  <  tfinal) 

disp( ‘SIN G UL A R ITY  LIK ELY. ’) 
t 

end 


end; 


%Display  the  venting  duration  and  total  mass  vented  on  the  screen 

disp(‘Vent  time(sec)  Total  vented  massfein)’) 

[Vtm  Tmvtd] 


%form  strings  for  placement  on  plot 


stl— ‘Duration  of  Venting  =’  ;st2=mnn2.str(Vtin);  st3=’  sec.’; 
st4=‘Total  Vented  Mass  =’;st.5=num2.str(Tmvtd);  st6=’  gins’; 
strl=[stl  st2  st:3];str2=[.st4  st5  st6]; 
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7oPLOT  THE  OUTPUT  TIME  HISTORIES 

%Note:  The  gtext(‘  ’)  items  are  for  placement  with  the  mouse 


if  kplt(l,l)  >  0, 

plotyy(toTit,oxout, tout, 0 lout, ‘Overpressure  Impulse  Ratio’);grid; 

title(‘EXAMPLE  FIG  1:  BLAST  VENTING-  OVERPRESSURE  AND  IMPULSE  HISTORIES’) 

xlabel(’Tinie,  sec’) 

ylabel( ’Overpressure  Ratio’) 

gtext(’Overpressure’); 

gtext(’Ovp  Impulse’); 

pause; 

figure 

else, 

disp(’No  Overpressure  or  Overpressure  Impulse  plots’) 
end; 


if  kplt(l,2)  >  0, 

plotyy(tout,DRout,tout,TRout, ‘Temperature  Ratio’);grid; 

title(‘EXAMPLE  FIG  2:  DENSITY  RATIO  AND  TEMPERATURE  RATIO  HISTORIES’) 

xlabel(’Time,  sec’) 

ylabel(’Density  Ratio’) 

gtext(’Deusity’); 

gtext( ’Temperature’); 

pause; 

figure 

else, 

disp(’No  Density  Ratio  or  Tem])erature  Ratio  plots’) 
end; 


if  kplt(l,3)  >  0, 

plotyy(tout,VELout,tout,MEout,‘Mach  Number  of  Efllux’);grid; 

title(‘EXAMPLE  FIG  3:  EFFLUX  VELOCITY  AND  MACH  NUMBER  HISTORIES  ’) 

xlabel(‘Time,  sec’) 

ylabel(‘Efflux  Velocity,  cm /sec’) 

gtext(‘Velocity’); 

gtext(‘Mach  No’); 

pause; 

figure 

else, 

disp(‘No  Efflux  Velocity  or  Mach  Number  plots’) 
end; 


if  kplt(l,4)  >  0, 

plotyy(tout,MSDout,tout,TOTMout,‘Total  Efflux  Mass,  gms’);grid; 

title(‘EXAMPLE  FIG  4:  EFFLUX  MASS  FLOW  RATE  AND  TOTAL  MASS  HISTORIES’) 
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xlabel(‘Time,  sec’) 

ylabel(‘Mass  flow  rate,  gin/sec’) 

gtext(‘Mass  rate’); 

gtext(*Total  Mass’); 

gtext(strl);%Prints  Venting  Duration 

gtext(str2);%Prints  Total  Mass  vented 

pause; 

figure 

else, 

disp(‘No  Mass  Flow  Rate  or  Total  Efflux  Mass  plots’) 
end; 


if  kplt(l,5)  >  0, 

plotyy(tout,DVout,tout,-DTout, ’Relative  Heat  Loss’);grid; 

title(  EXAMPLE  FIG  5:  RELATIVE  VOLUME  CHANGE  AND  HEAT  LOSS  HISTORIES’)’ 
xlabel(‘Tiine,  sec’)  '' 

ylabel(‘Relative  Volume  Change’) 

gtext( ‘Volume  Change’) 

gtext(‘Heat  Loss’) 

pause; 

else, 

disp(‘No  Relative  Volume  Change  or  Relative  Heat  Loss  plots’) 
end; 


function  xp=intair8(t,x,V,xxO,ssmax,area,TR,TRcut,gam) 


gm=gam;gl=gm-1.0;g2=0.rj*(gm+1.0)/gm;gm3=0.5*gl/gm; 

g3r=0.5/g:3;g4=-1.0/gm;gchok=(gm*g2)-g3r; 

gchl=((g2+gm)"(-g2*g3r))*((0.5*gl)-0.5); 

%  _ _ _ _ 

%It  IS  assumed  that  initially,  the  pressure  rise  due  to  explosion  in 
%the  chamber  is  instantaneous; 

% _ 

%  Choked  flow 


if (x/gchok)  >  1, 
AR=area; 


%No  choking 
elseif  (x/gchok)  <1, 

AR=area*((g2*gm)*(g2*gm/gl))*((2/gl)-0.5)*(x*(-g2))*((x-  (2*g3)-1.0)'0.5); 

else, 
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AR=area; 

end; 

xssmx=(xx0‘(-g3))*ssmax; 


%  Time  derivative  of  chamber  pressure 


if  TR>TRcut 

xp=-gm*((g2*gm)*(-g2*gin/gl))*(xssinx*(AR/V))*(x‘(l+g3)); 

elseif  TR<TRc.ut, 

xp=0; 

else, 


xp=0; 

end; 
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Figure  B-l.  Blast  venting  —  overpressure  and  overpressure  impulse  histories 


Figure  B-2.  Density  ratio  and  temperature  ratio  histories. 
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Overpressure  Impulse  Ratio 


Figure  B-3.  Efflux  velocity  and  mach  number  histories. 
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Figure  B-4.  Mass  flow  rate  and  total  efflux  mass  histories. 
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APPENDIX  C 
LIST  OF  SYMBOLS 


Roman  symbols 

a  =  Thermal  difFusivity  of  wall  material 

Co  =  Vent  cross-sectional  area 

ARq  —  Initial  area  of  chamber  wall  surface 

B  =  A  constant  in  computing  volume  change  due  to  thermal  expansion 

Cuj  =  Compressional  wave  speed  in  wall  material 

coag  =  A  reference  sound  speed 

cog  =  Initial  (maximum)  sound  speed  in  chamber 

cocr  =  Initial  critical  sound  speed 

C  =  Specific  heat  of  the  wall  material  of  chamber 

Cp  =  Specific  heat  of  air  at  constant  pressure 

/  =  Fraction  of  simulant  fluidized 

liiQ  =  Initial  (maximum)  Mass  Rate  of  Efflux 

Mq  =  Initial  Mach  Number  of  Efflux 

Mex  =  Mass  of  explosive 

Mh  m  —  Mass  of  loose  hazardous  material  simulant  inside  the  chamber 
Q  =  Gain  (or  loss)  of  heat  energy 
Q  =  Heat  flow  rate 

=  Increment  in  gain  or  loss  of  heat  energy 
po  =  Initial  (maximum)  pressure  of  chamber  gases 
POa  =  Reference  (ambient)  pressure  (=  1.01325*10®  dynes/sq.  cm) 

Pgit)  =  Quasistatic  pressure  in  the  chamber  at  time  t 
i2o  =  Initial  (maximum)  density  i-atio  (=  PogIPoa) 
t  =  Time  after  explosion 
6t  =  Time  increment 

To  =  Initial  (maximum)  temperature  of  chamber  gases 

Toa  =  Reference  (ambient)  temperature  (a  constant,  e.g.  298.15°  K) 

~  Quasi-static  temperature  of  chamber  gases  at  time  t 
u  =  Convection  coefflcient  for  wall  material  and  chamber  g8ises 
Uq  —  Initial  (maximum)  efflux  velocity  through  vent 
K)  =  Initial  volume  of  chamber 
w  =  Impedance  of  wall  material 
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W  =  Yield  of  explosion 

xog  =  Initial  (maximum)  pressure  ratio  in  chamber 

^cg  =  Cutoff  for  pressure  ratio 

Xg{t)  =  Normalized  queisi-static  pressure 

Greek  symbols 

Po  =  Constant  in  heat  loss  by  conduction 
X  =  Choking  pressure  ratio 

€  =  Coefficient  of  linear  expansion  of  wall  material 
7  =  Specific  heat  ratio  for  chamber  gases  (a  constant,  e.g.  1.35) 

K  =  Thermal  Conductivity  of  wall  material  of  chamber 
A  =  Constant  in  volume  change  adjustment  to  pressure 
POa  =  Ambient  air  density 

POc  =  Initial  density  of  chamber  contents  without  simulant 
POcr  —  Critical  density 

Pog  =  Initial  density  of  chamber  gases  including  fluidized  simulant 

Pg(t)  =  Density  of  chamber  gases  at  time  t 

Prv  =  Density  of  wall  material  of  chamber 

00  =  Initial  (maximum)  temperature  ratio  in  chamber 

©cp  =  Cutoff  for  temperature  ratio 

Qg{t)  =  Temperature  ratio 

Qgi  {t)  =  Interim  temperature  ratio  ^  =  Film  coefficient 
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